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Abstract. 

We study a general mass transport model on an arbitrary graph consisting of L 
nodes each carrying a continuous mass. The graph also has a set of directed links 
between pairs of nodes through which a stochastic portion of mass, chosen from a 
site-dependent distribution, is transported between the nodes at each time step. The 
dynamics conserves the total mass and the system eventually reaches a steady state. 
This general model includes as special cases various previously studied models such 
as the Zero-range process and the Asymmetric random average process. We derive a 
general condition on the stochastic mass transport rules, valid for arbitrary graph and 
for both parallel and random sequential dynamics, that is sufficient to guarantee that 
the steady state is factorisable. We demonstrate how this condition can be achieved 
in several examples. We show that our generalized result contains as a special case 
the recent results derived by Greenblatt and Lebowitz for c?-dimensional hypercubic 
lattices with random sequential dynamics. 
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1. Introduction 

Diverse physical phenomena such as traffic flow P, force propagation in granular 
media [21 E], clustering of buses |5, aggregation and fragmentation of clusters [S], 
phase separation dynamics [Hj, shaken granular gases [3 Ej and sandpile dynamics jH] 
share one common feature-their microscopic dynamics involves stochastic transport of 
'mass', or some conserved quantity, from one point in space to another. To simplify 
analysis, continuous space is typically replaced by (or "binned" into) discrete sites. 
Several such lattice models of stochastic mass transport have been introduced and 
studied, most notably the Zero- Range Process (ZRP) jUIl IHl [T2] , and the Asymmetric 
Random Average Process (ARAP) ^3 El- These models are deflned by specifying 
the microscopic dynamics, i.e. the basic stochastic rules for mass transport. Given the 
dynamics, there are two principal theoretical issues: (i) to identify the steady state if 
there is any, i.e. to flnd the invariant measure and (ii) once the steady state is known, 
to understand various physical properties in the steady state, e.g. the phenomenon of 
'condensation' that happens when a flnite fraction of the total mass condenses onto a 
single site [T5] . 

It turns out that the step (i) itself is often very difficult and the exact steady state 
is known in only very few cases JT]. In many of these known cases, the steady state is 
factorisable. This means that the steady state probability P({mj}) of flnding the system 
with mass mi at site 1, mass m2 at site 2 etc is given by a product of (scalar) factors 
f{mi) — one factor for each of the L sites of the system — e.g. for a homogeneous 
system where all sites i have equivalent connectivities 



where Z{M, L) is a normalisation which ensures that the integral of the probability 
distribution over all conflgurations containing total mass M is unity, hence 



Here, the 5-function has been introduced to guarantee that we only include those 
configurations containing mass M in the integral. The single-site weights, f{m) are 
determined by the details of the mass transfer rules and for a heterogeneous system 
may depend on the site i. 

The advantage of having a factorisable steady state is that the step (ii) mentioned 
above is often easier to carry out explicitly. This has been demonstrated recently by an 
exact analysis of the condensation phenomenon that occurs in a general class of mass 
transport models [T^ll6j . This raises a natural question: when does the steady state in 
these mass transport models factorise? This issue was recently addressed in the context 
of a sufficiently general 'mass transport model' in one dimension, that includes, as special 
cases, the previously studied ZRP, ARAP and the chipping model [HI- In this model 
a mass rui resides at each site i oi a, one dimensional lattice with periodic boundary 




(1) 
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conditions. At each time, a stochastic portion, rhi < of the mass mi at site i, chosen 
from a distribution (y9(mj|mj), is chipped off to site i + 1. The distribution ip{rh\rn) 
was called the chipping kernel and by choosing its form appropriately one can recover 
the ZRP, the ARAP and the chipping model of 5J. Even though the model above is 
defined in discrete time where all sites are updated in parallel, by appropriately choosing 
the chipping kernel it is easy to study the continuous time limit, which corresponds to 
a random sequential update sequence, as a special case [T7j. Similarly, one can also 
recover, as a special case, the model with only discrete masses as in ZRP. Thus the 
discrete time dynamics generalises continuous time dynamics but a continuous mass 
variable generalises discrete mass. 

A natural question, first addressed in ^7j, is what should be form of the chipping 
kernel ip{fh\rn) for the the final steady state to be factorisable. In that study, it was 
proved that the necessary and sufficient condition for a factorised steady state in the 
one dimensional directed case defined above is that the chipping kernel is of the form 

/ ~ I X uim) vim — m) 

V{m\m) = \ '\ ' (3) 
Jq dm u\m)v[m — m) 

where u{z) and v{z) are arbitrary non-negative functions. Then the single-site weight 
in Eq. (|T]) is given by 

I'm 

f{m)= I dmu{m) v{m — m). (4) 



Furthermore, given a chipping kernel (f{m\m), sometimes it is hard to verify explicitly 
that it is of the form Q and thereby to identify the functions u{m) and v{m) in order 
to construct the weight f{m) in Eq. (^. This problem was circumvented by devising a 
test I8j to check if a given explicit ip{m\m) satisfies the condition Q or not. Further, if 
it "passes this test," the weight /(m) can be found explicitly by a simple quadrature [IHj. 
Finally, for any desired function f{m), one can construct dynamical rules (i.e., ip) that 
will yield /(m) in a factorised steady state. 

It was further demonstrated in Ref. JZj that the corresponding necessary and 
sufficient condition in the case of random sequential dynamics in continuous time can be 
easily obtained, by taking a suitable limit, from the condition for the parallel dynamics 
manifest in Eq. Q. This is done by choosing the chipping kernel as 



(p{m\m) 



1 — dt / a(m\m) dm 







6{m) + a{m\m) dt (5) 



for small time increment dt and 6{z) is the Dirac delta function. The function a{m\m) 
denotes the 'rate' at which a mass m is transferred from a site with mass m to its right 
neighbour. Note that the form in Eq. (jS)) ensures the normalization, (f{m\m) dm = 1. 
Then, the necessary and the sufficient condition for factorisable steady state, derived 
from the more general condition in Eq. (jS)), is that the rate a{m\m) must be of the 
following form ^\ 
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where x{z) and v{z) are arbitrary non-negative functions. The corresponding steady 
state weight is then simply, f{m) = v{m). 

The condition (jH]) for factorisabihty in the mass transport model was derived only 
in one dimension and also only for unidirectional mass transport (from site i to site 
i + 1). A natural question, therefore, is whether one can generalise this condition to 
higher dimensional lattices, or to arbitrary graphs where mass transport can take place, 
in general, between any pair of sites i and j. Recently, Greenblatt and Lebowitz were 
able to derive a sufficiency condition ^Hj for factorisabihty in the mass transport model 
with nearest neighbour mass transport on a regular d-dimensional lattice with periodic 
boundary conditions, but considered only the case of random sequential dynamics. They 
showed that if a;q(m|m) dt is the probability of mass m < m being chipped off a site with 
mass m to a nearest neighbour in the direction q (there being 2"^ nearest neighbours on 
a hypercubic lattice in d dimensions), then the sufficient condition for factorisabihty is 
a direct generalization of the condition in Eq. ®, namely that the rate function must 
be of the form 



for each q, where Xq{z) for each q and v{z) are arbitrary functions. The steady state 
weight is simple, /(m) = v{m). However, it was not possible to prove that the condition 
((Zj) is also necessary [TH], except in the case of generalized Zero-range processes. 

The purpose of this paper is to derive a more general sufficiency condition, valid for 
arbitrary graphs where the mass transport takes place not necessarily between nearest 
neighbours and for the more general case of parallel dynamics. Our results boil down 
to equations (j^^ and pip. The former yields a sufficiency condition and the latter an 
additional consistency condition which must be satisfied. 

In the special case of random sequential dynamics on a regular hypercubic lattice 
with nearest neighbour mass transport, our sufficiency condition reduces to the one 
derived by Greenblatt and Lebowitz. Our results, however, are considerably more 
general. 

The paper is organized as follows. In Section 2, we define the mass transport model 
on an arbitrary graph. In Section 3, we derive a sufficient condition on the chipping 
kernels, valid for parallel dynamics on an arbitrary graph Q, that would gaurantee that 
the steady state is factorisable. We show that there are some additional consistency 
conditions that need to be satisfied in general and we demonstrate explicitly how these 
conditions are satisfied in several examples. In Section 4, we extend our approach to 
random sequential dynamics on an arbitrary graph. We conclude with a summary and 
discussion in Section 5. 

2. The Mass Transport Model on an Arbitrary Graph 

We consider a fixed arbitrary graph Q consisting of L nodes labelled i = 1, 2, 3, ■ ■ ■ , L 
and a set of directed links or channels between certain pairs of nodes. At a given time 




(7) 
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t, a node i has mass mj > where mi are continuous variables. We consider discrete 
time dynamics where at each step the masses at all the nodes are updated in parallel 
according to the following rules. We first define a (L x L) mass-transfer matrix [/i] as 
follows. An element fiij of the matrix [/i] is identically zero at all times if there is no 
directed link from site i to site j on Q. If there is a directed link from i to j, then 
Hij > is a non-negative stochastic variable that represents the mass transferred from 
site i to site j during one update. In Fig. 1 we give an example which we refer to for 
illustrative purposes throughout this section. The diagonal element fiu represents the 
mass that stays at site i at the end of the single update. We assume that the dynamics 
of mass transport conserves the total mass. Thus if {mi,m2, . . . ,rnL} represents the 
masses before the update, by virtue of mass conservation, the row sums of the matrix 
[fi] are given by, J2j t^ij = (see Fig. 1). Similarly, the column sum corresponding to 
a node i, J2j f^ji = ""^i represents the mass at site i just after the update. 

Note that, by definition, some elements of the matrix [/x] are permanently zero 
(when there is no directed channel available for mass transfer between a pair of sites). 
On the other hand, when there is an available channel from itoj, the matrix element /i^ 
is a stochastic variable which is chosen in the following way. For each node i, we define 
a generalized 'chipping kernel' ipi ({/ijj}|mj) that represents the joint distribution of the 
masses transported from site z in a single update. Here the set {fiij} runs over only 
those sites j which are connected to i via a directed link, and in addition it includes the 
diagonal element fin. In other words, {fiij} is simply the set of non-zero elements in the 
i-th row of the mass-transfer matrix [/i]. For example, in Fig. 1 where we have a graph 
of four nodes along with the directed links, we need to define four chipping kernels as 
follows: ipi {fill, flu, fiulmi), (/^2i, /i22, yU23|"^2), V's (Ai3i,/i33|"^3) Siiad ip^ {fi43, fiulm^) 
respectively. The chipping kernels must be normalized to unity at all sites i, i.e. 



where it is implied that the index j in the sum as well as in the product runs over 
all the sites that i feeds into, including i itself. Thus in this model the amount of 
mass transported from a given site i in one update depends only on that site i, and 
not, e.g. on the destination sites to which the mass is transported. Also, we assume 
that the chipping kernels ipi do not contain the time t explicitly. Note also 

that this chipping kernel ipi {{fiij}\mi) generalises that of (TTHIH] sufficiently to include 
hypercubic lattices ^H] or more complicated graphs. 

The chipping kernels thus specify the dynamics, i.e. the mass update rules. Given 
these kernels, we next ask what is the steady state joint probability distribution of masses 
m = {mi, 7712, ■ ■ ■ , mi}, i.e. P [m, t oo). In particular, our goal is to determine the 
properties of the chipping kernels ipi {{fiij}\mi) required in order to guarantee that the 
steady state joint probability distribution is factorisable, i.e. of the form 
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Figure 1. An example graph with four nodes labeUed 1, 2, 3 and 4 and directed 
links between certain pairs of nodes. The associated (4 x 4) mass-transfer matrix [/i] 
is shown, whose element /i^ denotes the stochastic mass transferred from site i to 
site i in one single update, provided there is a directed link between the two sites. 
If there is no directed link, the corresponding matrix element is always identically 
zero. The diagonal element [in is the amount of mass that stays at site i during the 
update. The row sum and the column sum associated with any node i, /x^ = rrii 
and "Ylij P'ji — '^iJ represent respectively the mass at i before and after the update. 



where the normahzation constant is given by 
Z(M,L) = n 

i=i 

Besides, if the steady state factorises as in Eq. 0, we would also like to know the single- 
site weights fi{mi) in terms of the prescribed chipping kernels (pi ({/ijj}|mj). Note that 
on an arbitrary graph Q, the single-site weights fi{mi) are, in general, different from site 
to site. Hence there is an additional subscript i in fiinii). On a homogeneous graph, 
where all sites have equivalent connectivities, the weight function /(m) does not depend 
on the site i explicitly as in Eq. ((H). 

3. Factorisable Steady State on an Arbitrary Graph 

Since the dynamics conserves the total mass, at any time t we can write P (m, t) oc 
F{rn,t)5 (Yli^nii — M) where M is the total mass and F{m,t) is the unnormalized 
weight at time t. Below, we will first write down the general evolution equation of 
the weight F{m,t) under the mass transport rules prescribed by the chipping kernels. 
While the notations that we will use for a general graph Q may seem a bit complicated, 
it is instructive to keep the simple example in Fig. 1 in mind and use it as a guide to 
the general notations. 



dniifi 



ITLi 



M 



(10) 
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Let us consider a single update from time t to time t+ 1. Let m = {mi,m2, ■ ■ ■ , mi} 
denote the masses at time t before the update and rn/_ = {m'^, m'g, ■ ■ ■ , m'j^} denote the 
masses at time t + 1 after the update. In terms of the elements of the mass-transfer 
matrix [/i], we thus have, mj = ^^./ijfc and m'- = Ylikl^ki for all i. The master equation 
for the evolution of the weight then reads 

F{m[_,t+l)= ( 

L 



n 

2=1 



(11) 
F(m, t) 



k k 

where the product over j runs over only the sites to which the site i feeds into, i.e., when 
there is a directed link between sites i and j (note that this set includes the site i itself). 
For other sites that are not connected to i by a directed link, the corresponding matrix 
element = identically (see Fig. 1) and hence they are not integration variables in 
Eq. (HH). 

Our strategy now is to assume that the steady state factorises and to determine a 
sufficent condition for this assumption to hold. First, we take the t — oo limit on both 
sides of Eq. f|TT|l and assume that the steady state weight factorises 

F{rn) = l[f,{m,) (12) 



and write 



thus 



fi{mi)ipi {{^ij}\mi) = Vi{{fXij})5{mi - ^yUy) 



(13) 



n/(^^)=n 

j=i 1=1 

Now we Laplace transform this equation to obtain 

L L 

n^?(^)=n 



i=l i=l 

where we have defined gi{si) = fi 



/•oo /• 

/ druiYl di2ijVi{{i2ij})6{mi -^fHk) e""'^^'^*' 
Jo , J J, 



.(15) 



m)e' 



'Si m 



dm. Next we trivially perform the 



integration over the mj variables on the rhs of Eq. ()15j) and rearrange the Si to get 



n^(^)=n 



i=l 



i=l 



Y[ j dnijVi{{iJij})e ^ 



llX,{{s,h). (16) 



where {sk}i indicates that k runs over the sites to which i is connected by a directed 
link, including the site i itself. For example, for the graph in Fig. 1, we have 

giisi)g2{s2)g3{s3)g4iSi) = Xi{si, S2, S4) X2{si, S2, S3) X3(si, S3) X4(s3, S4).(17) 

One solution (not the most general as we shall discuss in section El ) to ()16|) is 

X.{{s,h) = llK,,is,) (18) 
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(19) 



As before, the product in Eq. (fTH|l runs over only those sites j that the site i feeds into 
(including itself), whereas in Eq. (fT^ the product runs over only those sites j that feed 
onto i (including the site i). For example, for the graph in Fig. 1, we will have 

Xi{si, S2, S4) = Ku{si) Ki2{s2) -^'14(54) 

X2(si, S2, S3) = K2l{Si) K22{S2) -?^23('S3) 

X^isi, S3) = i^3i(si) ^^33(53) 

^4(53,54) = 7^43(53) i^44(s4) (20) 

and correspondingly 

gi{si) = Ku{si) K2iisi) i^3i(si) 

92(82) = K 12(3 2) K 22(3 2) 

5(3(53) = ^^23(53) ^33(53) ^43(53) 

54(54) = K14{S4)K44^{S4). (21) 

Since the rhs of Eq. (fT^ is a product of Laplace transforms, its inverse Laplace 
transform fi{mi) therefore has the convolution form 



fiim) 



n 



*3 



(22) 



where the notation indicates a product over sites j feeding into i which is a multiple 
convolution and the function Vji{m) is the inverse Laplace transform of Kji{si), i.e. 



'da 



(23) 



Again, going back to Fig. 1, we will have, for example, by inverting Eq. (pT|l the 
following convolution for the site labelled 1 



fiimi) = j vu{mi - (T21 - 0-31) ^^21(0-21) ^^3i(cr3i)'icr2i da^i 



(24) 



and similarly for the sites labelled 2, 3 and 4. 

Having obtained the form of fi{m), let us ask what does this imply for the chipping 
kernels. To see that, we go back to the definitions of Xi in Eq. p6|) . Substituting the 
ansatz for Xi in Eq. (|T8|) on the rhs of Eq. |T6|) we see that for each i 



n / dfi.,V.{{fi.,})e-^^^''^^''=llK,,{s,) 

For example, for the site 1 in Fig. 1, we have 

j d/iii d^i2 dfJ^u ^i(/iii, /U12, /iu) e~'^''''' 



'S2 M12— S4 fJ-14 



(25) 



K,,{s,)K,2is2)K,4{s4),i26) 
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and similarly for the sites 2, 3 and 4. Eq. immediately implies for the general case 

'P^i{^^^l})=lly^Ml) (27) 

j 

where j runs over the sites which site i feeds into and Vij is the Laplace inversion of Kij 
as in Eq. (j2Sl)- For example, for the site 1 in Fig. 1 we have 

-pL(/iii,/ii2,/ii4) = Vn{iJ,n) VuifJ^u) VuilJ^u)- (28) 
Substituting Eq. (jTfjl in Eq. (fT!?jl we get the following expression for the chipping kernel 

II N Ylj'^ijif'ij) 

•"'""""'"■'^ in.. ^"..] ('».)■ 

In the numerator the product is over sites j which site i feeds into via a directed link 
connecting i to j (this includes site i itself), and one has the constraint J2j f^ij — 
In the denominator, on the other hand, the product is over sites k feeding into site i 
(including once again the site i itself). Going back again to the example in Fig. 1, we 
have for the site 1 the following chipping kernel 

/ I X ^Ii(/iii)^^i2(/ii2) ^14(^14) ,on^ 

<^l(/^ll,/il2,/il4P^l) = 7 7 ^ 7 7 7 ri 1 (30) 

j Vu{mi - 0-21 - 0-31) ^21(0-21) ^^3l(0^3l)"Cr21 dasi 

where fin + fiu + /ii4 = nii. Similarly one can easily express the chipping kernels for 
the other sites 2, 3 and 4 in terms of the functions Vij. 

However this is not quite the end of the story since the chipping kernel (pi{{fiij}\mi) 
must obey a key consistency condition, namely its normalization: Eq. (jSJ. This implies 
from Eq. that for each site i 



n 



(m) 



n 



(m) (31) 



Again, the convolution on the Ihs is over sites j feeding into site i but the convolution 
on the rhs is over sites j which i feeds into. For example, for the site 1 in Fig. 1, this 
consistency condition reads 

VwivTLi - /ii2 - /Um) Vyi{^\2) V\^{fi\^) dj^u dj^u. (32) 



One can write down similar consistency conditions for sites 2, 3 and 4 also. 

In summary, a sufficient condition for factorisation is that the chipping kernels 
are product form as in Eq. ()29j) where Vij are non-negative functions that satisfy the 
consistency conditions in Eq. ()31|). but otherwise are arbitrary. If this condition holds, 
then the steady state will factorise with single-site weight fi{m) given by the convolution 
formula in Eq. (j221)- These conditions form the central result of this paper. The 
conditions arise from making the ansatz in Eq. (fTHj) which immediately implies Eq. 
(j221) for fiimi) and Eq.(j2Z|) for Vi{{iJ,ij}). These two equations, when substituted in 
Eq. ()13p. gives Eq. ()29p for the chipping kernel. Finally, Eq. ()31|) arises from the 
requirement that the chipping kernel in Eq. is normalized to unity. 
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Condition ()29pHH) is certainly a sufficient condition. Whether it is also a necessary 
condition remains a non-trivial issue. In the special case of one dimensional graph with 
unidirectional transport, it was explicitly proved in Ref. that this is also a necessary 
condition. Similarly, this is true (see below) in the case of a complete graph where every 
site is connected to every other site. However, it is unlikely, in general, that ()29|l is 
also necessary. In section we provide an explicit example showing solutions to (fTBjl 
which are not of the form ()29p. Although we are unable, for a technical reason to be 
discussed later, to prove that this class of solutions are also "legitimate," we feel that 
valid solutions different from ()29|) should exist. Beyond this simple example, finding the 
most general valid solution to Eq. (|TH|) for the arbitrary graph seems to pose significant 
challenges. 

3.1. Specific Examples 

Let us now discuss several special cases where one can directly verify that the consistency 
conditions in Eq. (jSH) are satisfied automatically by the functions Vij due to the nature 
of the underlying graph Q, thereby guaranteeing factorisability provided the chipping 
kernels are of the product form as in Eq. ()29|) . 

Example 1: Mass transport in one dimension over extended range: Consider 
a one dimensional latice with periodic boundary conditions where mass from a given 
site can get transported up to a range / to the left or to the right. In this case 

{{fJ'ij}\m) ■ ■ (33) 

and we have a translationally invariant chipping kernel i.e. (pi{-\m) is the same function 
of (2/ + 1) arguments for each site. Then the consistency condition (jHlll reads 

which is automatically satisfied due to the translational invariance Vi^i+dif) = Vi_d,i{.l^)- 

Example 2: Symmetric chipping rule: It is easy to see that if the chipping rules are 
symmetric, i.e. fjj(/i) = Vji{fi), then the consistency condition holds automatically. 

Example 3: A more general chipping rule: A more general way to meet the 
condition ()3ip is to assume that the chipping rules are such that at each site i, for every 
directed link {i j) out of site i, there is an incoming directed link {k i) to i from 
some other site k with Vij{fi) = ffcj(/i). For example, in the 1-d example above, this is 
achieved through Vi^i+d{fJ') = Vi-d,i{ft). Similarly, one can achieve this on a hypercubic 
lattice in arbitrary dimensions and on any graph where the number of links out of a site 
equals the number of links into it. 

Example 4: Complete graph with permutationally invariant chipping kernel: 

So far we have proved that Eq. ()29|1 is a sufficient condition for factorisability on 
any arbitrary graph provided we can satisfy the consistency condition (jM]) for each 



Factorised Steady States in Mass Transport Models on an Arbitrary Graph 



11 



site. As discussed before, to prove that this condition is also necessary seems hard 
for a general graph. However, one can prove this for a complete graph where every 
site is connected to every other site via a directed link and the chipping kernel is 
permutationally invariant. This can be proved as follows. On a complete homogeneous 
graph the functions gi{s) = g{s) in Eq. (fTHjl do not depend explicitly on the site index 
i. Moreover, since every site is connected to every other sites and using the fact that 
the chipping kernel is permutationally invariant Eq. ()16p simply becomes 

9isi)g{s2) ■ ■ ■ g{sL) = [X{si, S2, ■ ■ ■ , sl)]^. (35) 

Clearly X{si, S2, ■ ■ ■ , sl) = K{si)K{s2) ■ ■ ■ K{sl) with K{s) arbitrary, is one solution 
to Eq. (jnni) and thereby g{s) = [K{s)]^. As in the general case, this will then lead to 
the sufficiency condition (j29|l . To prove that X{si, S2, ■ ■ ■ , sl) = K{si)K{s2) ■ ■ ■ K{sl) 
is also the most general form of the solution that one can write down for Eq. (jSSl), we 
take logarithm on both sides of Eq. and then take derivatives with respect to Sj 
and Sj with i j. This gives 

dHnX _^ 

dsidsj 

for any i ^ j. It is then easy to see that the most general solution of the partial 
differential equation is of the form, X(si, S2, ■ ■ ■ , s^) = Ki{si) K2{s2) ■ ■ ■ Kl{sl) 
where Ki{s) are arbitrary functions. Since the graph is homogeneous, we also 
have Ki{s) = K{s) independent on the site index i. Thus X{si, S2, ■ ■ ■ , sl) = 
K{si) K{s2) ■ ■ ■ K{sl) with K{s) being an arbitrary function, is the only solution of 
Eq. (jH^j) that respects homogeneity. Since this solution finally leads to the sufficiency 
condition, the uniqueness of its form guarantees then that Eq. is both necessary 
and sufficient. Note that the consistency condition ()31|) is automatically satisfied in this 
case. 



4. Random Sequential Dynamics 

The sufficiency condition in Eq. (j^^ and the associated consistency condition in 
Eq. (jnH), derived above for parallel update dynamics in discrete time, can be easily 
extended to the case of random sequential dynamics. This can be achieved by letting 
the probability of the chipping event in a time step oc dt so that, to leading order in dt for 
small dt, at most one chipping event can occur in the whole graph Q per update, i.e. the 
chipping events occur sequentially one per update. In addition, taking dt ^ one can 
obtain the continuous time limit where chipping events occur with 'rates' per unit time. 
Thus, the corresponding sufficiency condition for the random sequential dynamics will 
be specified in terms of the chipping 'rate' kernels, rather than the chipping probablity 
kernels (pi in parallel dynamics. To translate the sufficiency condition in Eq. ()29|1 valid 
for the probability kernels into one that is valid for 'rate' kernels, we first redefine the 
functions Vij{^ij), for all i 7^ j, in the following way 

VijifJ^ij) = + Xij{fiij)dt (37) 
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where Xij{fiij) are arbitrary functions. The diagonal functions Vulfiu) are left unchanged. 



With this re-definition of Vij, the steady state weight in Eq. 



(m) + dt 



^m(^ l^ji) Xjii^Uji) dfi 



'Jl 



becomes 



+ 0{dt^ 



(3^ 



where the sum over j in the second term runs over all sites j ^ i that feed into 
site i on Q. Using the re-defined Vij in Eq. ()37|) one can similarly rewrite the 
chipping kernels in Eq. ()29|1 . Since the diagonal elements play a special role, it is 
convenient to redefine the chipping kernel only in terms of non-diagonal elements, i.e., 
ipi {{nij}\mi) = ipi ({/ijj}'|mj) where {^ij}' denotes the set of matrix elements in the row 
i without the diagonal element ^u. We are allowed to get rid of the diagonal element 
using the row sum, Ylj f^ij ~ ■^'^^ example, for the graph in Fig. 1, we will rewrite, 
^lifJ'ii, fJ'12, fJ'ul^i) = fJ'i'il^i)- Substituting Eq. (pTTjl in Eq. (jSHI) and taking 



the limit dt 



^ 0, one gets 

{{i^ijyim 



dt 



Viiim. 



+ dt 



* Vii] {m, 



V. 



+ 



Vidnii 



ik 



+ 0{dt^) (39) 



where [x * y\{jn) = x{a)y{m — a)da denotes the convolution integral. The notation 
in Eq. ()39p may look a bit complicated, but actually it's rather simple. For example, 
for the graph in Fig. 1, the kernel in Eq. for the site labelled 1 reads. 



dt 



(40) 



(Ki * a;3i](mi) + [vu * a:2i]("^i)) 



-f- dt 



vii{mi) 

Xl2{tl'12)Vu{mi - flu) 



. , XiA{fllA)Vll{mi - flu) 

^(/Um) H 7 — 7 d[iii2) 



viiimi) viiimi) 

Taking dt ^ Q limit in Eq. ()39|) we then obtain the continuous time limit where a 
portion fiij gets chipped off the site i with mass mj, to another site j which is connected 
to z by a directed link from i to j on with a rate Xij{fiij)vii{mi — fiij)/vii{mi). Also, 
we notice that in the limit dt 0, the steady state weight fi{mi) is given simply from 
Eq. 



fiinii) = Viiinii). (41) 

However, as in the case of parallel update, we need the functions Xij to satisfy certain 
additional consistency conditions. This is because the chipping kernel ipi ({/ijj}'|mi), 
when integrated over its arguments {fiij}' must give unity. Integrating Eq. (jH^jl then 
gives the required consistency condition that must be satisfied for each node i, 

^ [xij * Vii] (mi) = ^ [xji * Vii] (mi), (42) 
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where the sum on the Ihs runs over all sites j that the site i feeds into (excluding the 
site i itself), but on the rhs the sum runs over all sites j that feed into site i (excluding 
the site i itself). As an example of this condition for the site labelled 1 in Fig. 1, we 
have 

[3^12 * vu]{mi) + [xu * vii]{mi) = [x2i * fii](mi) + [0:31 * Vu]{mi). (43) 

Similar consistency conditions can be wriiten down for sites 2, 3 and 4 in Fig. 1. 

Thus, in summary, for the mass transport model on an arbitrary graph Q with 
random sequential dynamics specified by the rates aij{fi\mi) of mass /i to be transferred 
from site i to site j (provided there is a directed link between i and j on ^), a sufficiency 
condition for the steady state to be factorisable with weights fi{mi) is that aij{fi\mi) is 
of the form 

"^^■^^'"^^^ = ^) ^''^ 

for all i, where Xij and vu are functions that must satisfy the consistency conditions in 
Eq. but otherwise arbitrary. Also, the corresponding steady state weight is then 

given simply by fi{mi) = Vii{mi). Note that though formally Eq. looks like a direct 
generalization of the one dimensional condition in Eq. the additional consistency 
conditions in Eq. P2|l are nontrivial to satisfy. For the one dimensional example in 
Eq. ©, the corresponding consistency condition, [xj^j+i * t>j](mj) = [xi-i^i * fj](mj) is 
automatically satisfied due to the translational invariance, Xj„i^i(/i) = Xi^i+i{fi). 

4-1- Specific example of hypercubic lattice 

It is easy to verify that the consistency condition in Eq. ()42|) is automatically satisfied in 
the four example cases of section 13.11 (as it should be since random sequential dynamics 
is just a limit of the discrete time case). Let us do this explicitly for the case where 
the graph ^ is a homogeneous hypercubic lattice with periodic boundary conditions and 
mass transfer takes place only between nearest neighbours. For a hypercubic lattice, 
from each site i there are 2'^ outgoing links to the 2'^ nearest neighbours of i. Similarly, 
there are exactly 2'^ incoming links to site i from its nearest neighbours. Then the 
sufficiency condition in Eq. ()44|1 can be written in a simplified notation 

ai,g(/i mi) = — (45) 

Vii{mi) 

where the index q = ±1 runs over the 2'^ directions {±ei, ±62, ■ ■ ■ , ie^}, ai^q{fi\mi) 
denotes the rate of transfer of mass fi from site i with mass mj in the direction q and 
Xi^q denotes the function associated with the link {i,i + q). Similarly, the the consistency 
condition in Eq. can be rewritten as 

^[xi,g * Vii]{mi) = ^ [xj-q * Vii]{mi). (46) 

q jgneighbours of i 

We next use the fact that the lattice is homogeneous, i.e. it is translationally 
invariant in all directions. Clearly then Xi-q^q{z) = Xi^i+q = gq{z) due to translational 
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invariance in the g-th direction. In that case the condition in Eq. is clearly satisfied 
at all i. Also, due to the translational invariance, the rate function aj^g(/i|mj) = aq{fi\mi) 
does not depend on the site index i. By the same requirement, vuimi) = v{mi). Thus, 
the sufficiency condition in Eq. ()44|) simply reads, 

aMm) = (47) 

where gq{z) and v{z) are arbitrary non-negative functions. The consistency conditions 
are automatically satisfied as proved above. The steady state single site weight /(m,) is 
simply, f{m.i) = v{mi), and naturally it does not depend on the site index i explicitly. 
The condition in Eq. (jTfj) is precisely that derived by Greenblatt and Lebowitz. Thus, 
our general sufficiency condition, valid for an aribtrary graph Q, recovers this special 
case when ^ is a homogeneous hypercubic lattice with nearest neighbour mass transport. 



4-2. ZRP on arbitrary graph 

We can also check that the steady state for the continuous time zero-range process on 
an arbitrary graph is recovered. The zero-range process involves discrete masses and is 
specified by the rates for a a unit of mass to hop from site i to j. In the case where 
these rates are given by 

aij{l\mi) = yijWi{mi) , (48) 

where Wi{mi) is the total rate for the unit mass leaving site i and yij is the probability 
that the random destination for a hop from site i is j, the steady state factorises with 
single-site weight 

Mm,) = (49) 

where Pi = Ylj^ti VjiPj the steady state probability of a single random walker moving 
from site i to j with probability i/ij [IIIIEI- 

To make the connection between the forms ()48|) and ()44|) we identify 



Xij[l) = yijPi and Wi{mi) = — - — . (50) 

Vii{rrii ) 

Inverting the latter equality to express vuimi) in terms of Wi yields, by virtue of (jUj), 
the single-site weight fl49|) . Finally the consistency condition ()42j) becomes 

'^Xij{l)vii{m - 1) = ^Xji(l)fii(m - 1), (51) 

and is satisfied with Xjj(l) = yijPi since we have 

^x.,(l) = ^a:,.(l). (52) 
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5. Conclusion 

In this work we have derived the sufficent condition ()29|1 along with a consistency 
condition (jHT|) for factorisation of the general (discrete time, continuous mass) mass 
transport model on an arbitrary graph. In this case the single-site weight is given 
by (j22I)- We gave in section ItTTl specific examples of geometries where the additional 
consistency condition associated with the sufficient condition is automatically fulfilled. 
Moreover on a complete graph with permutationally invariant chipping functions we 
showed that the sufficent condition is, in fact, also necessary. 

Of course a significant improvement would be to generalize condition ()29|31|) to a 
necessary and sufficient condition. To illustrate the challenges involved in accomplishing 
this task, we offer another simple example where we can derive a condition that is both 
necessary and sufficient. 

This example, a seemingly trivial generalization of the one in ^Tj E], involves a 
one-dimensional lattice with chipping to both nearest neighbour sites and we also assume 
that the chipping kernel is translationally invariant (/ = 1 in Example 1 of section ITTj) . 
Then, due to the translational invariance, (|T6|) becomes 

i i 

Using a similar procedure to that employed in Example 4 of section 13.11 — taking the 
logarithm of (j53|) then successive derivatives with respect to Si and sj — we find that 
the general solution to (j^Hj) is 

X{si-i,Si,Si+i) = Kii^i{si^i)Kii{si)Kii+i{si+i) ^}^''^' . (54) 

H{Si, Si+i) 

where Ku^i, Ku, Ku^i and H are arbitrary functions, independent of i. Inverting 
the Laplace transforms to give V^fiu-i, fin, ^u+i) would yield a considerably more 
complicated form for the chipping kernel ip^^a-i, fin, fiu+i) than (jSHI)- In principle, 
therefore, there could be a whole family of chipping kernels, generated by the choice of 
H, which give rise to the same steady state (i.e. the same single-site weight f{m)) as 
that for H = 1. Thus, H could be thought of as a "gauge function". 

In addition, for each choice of H, one has to ensure the consistency condition, 
namely that 

(p[fiii_i, fin, fiii+i) = — — r d[mi - [fiii-i + flu + fiii+i)) (55) 

is correctly normalized to unity. On the face of it, it may appear that this imposes a 
formidable constraint on the choice of H. However, we now show that this consistency 
condition does not impose any additional constraint on H. In other words, if the 
consistency condition is ensured for H = 1, then it is automatically satisfied for all 
other choices of H. To see this, consider the expression of the chipping kernel (j55|l . 
Now, the denominator f{mi) is, of course, independent of the choice of H. So, to prove 
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that the consistency condition of normahzation of (p does not impose any additional 
constraint on H, one has to prove that the integral 



/ (m) = J dfiii-idfiiidfiii+iV{fiii-i, fin, Hu+ijd (m - Hu^i - fin - fin+i) (56) 

is independent of the choice of H. Now, taking the Laplace transform of Eq. ()56|) one 
obtains 



where we have used definition (|TT))l in going from the first line to the second and Eq. 
()54|) in going from the second line to the third. But this expression on the rhs does not 
contain the "gauge function" H thereby proving that the integral of V is independent of 
the choice of H. Thus the consistency condition is automatically fulfilled for arbitrary 
choices of H as long as it is ensured for H = 1, which is indeed the case as shown in 
Sec, mil Example 1. 

A more serious constraint on H is that the inverse Laplace transform of ()54j) must be 
non-negative. The case H = 1 in (jH^ obviously imposes the trivial constraint that the 
inverse Laplace transform of the K^s be non-negative. It remains to be shown whether 
there is also a class of if 7^ 1 which lead to valid chipping kernels. If it could be shown 
that this class is non-empty, the one-dimensional / = 1 example we have discussed would 
show explicitly that it is not necessary for ip to be of the form ()29|). Mapping the extent 
of this class of H (if indeed it is non-empty) remains a daunting task. 

For an arbitrary graph, we have made some in-roads with similar considerations. 
Unfortunately, the most general solution to (fTT)|l remains difficult to formulate. 

Finally we note that having a factorised steady state opens the door for the study of 
condensation as in ^SlGSj- Thus one should be able to analyse condensation in various 
geometries or even on scale-free networks [22] ■ 
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